Topological phases and phase transitions on the square-octagon lattice 
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We theoretically investigate a tight binding model of fermions hopping on the square-octagon 
lattice which consists of a square lattice with plaquette corners themselves decorated by squares. 
Upon the inclusion of second neighbor spin-orbit coupling or non-Abelian gauge fields, time-reversal 
symmetric topological Z2 band insulators are realized. Additional insulating and gapless phases are 
also realized via the non-Abelian gauge fields. Some of the phase transitions involve topological 
changes to the Fermi surface. The stability of the topological phases to various symmetry breaking 
terms is investiaged via the entanglement spectrum. Our results enlarge the number of known 
exactly solvable models of Z2 band insulators, and are potentially relevant to the realization and 
identification of topological phases in both the solid state and cold atomic gases. 

PACS numbers: 71.10.Fd,71.10.Pm,73.20.-r 



I. INTRODUCTION 

Topological phases of matter have recently been the 
focus of intense theoretical and experimental effort. No- 
table amon g th em are the integer and fractional quantum 
Hall liquidSjlSS which cannot be understood in terms of 
the traditional description of phases based on symmetry 
breaking and local order parameters. While the quan- 
tum Hall states experimentally arise under strong mag- 
netic fields, a new paradigm of topological phased^ has 
emerged: the so-called topological band insulators (TBI) 
which occur in the presence of time-reversal symmetry 
and spin-orbit coupling.®® 

The TBIs are similar to the quantum Hall states in 
that they possess a gapped bulk spectrum and gapless 
edge states. However, a remarkable feature of TBI is 
that they can occur in both twcPHSl and three spatial 
dimension4i2ttl2l (while quantum Hall states are restricted 
to two dimensions). Following the initial predictions 
of a two dimensional TBI in HgTc quantum wells,'^ 
experi ment revealed this intriguing quantum phase of 
matter .EES! Not long afterwards, predictions were made 
for three dimensional compounda^^"^ and a number of 
these have now been verified experimentally.'^SHSs' Thus, 
TBI are now an experimentally established quantum 
state of matter. 

The classificatiorPSl of TBI is based on a Z2 
invariani^ * 10 12 27. j-^j. invariants, depending on the spa- 
tial dimension), rather than the Chern number used in 
the integer quantum Hall systems. The Z2 number turns 
out to be related to the parity of the number of gap- 
less edge modes (Dirac nodes) appearing on the surface 
of an insulator: Any odd number is topologically non- 
trivial and classifies the insulator as a TBI, while any 
even number can be shown to be adiabatically connected 
to the case with no gapless edge modes. For a band insu- 
lator, the Z2 n umber can be directly computed from the 
band structure.'^EEiESllZI Wjth a straightforward proce- 
dure in hand to classify band insulators as topological or 
"trivial" , the search is on to determine which models, and 
therefore which physical systems, are expected to reveal 



TBI physics. 

Clearly, real systems involve electron-electron interac- 
tions, but in many cases these can be treated a ccurately 
within a mean-field approximatiorP^ElHllEMg] j^^ which 

case the relevant physics comes down to single-particle 
band physics. Therefore, it is important to understand 
which features of a system lead to topological proper- 
ties in the band structure, whether that band structure 
is derived from a non-interacting model or results from 
the self-consistent treatment of an interacting problem. 
We emphasize that even in the absence of "microscopic" 
spin-orbit coupling, TBI can result in interacting sys- 
tems at the mean-field level by spontaneously generated 
spin-orbit coupling.l^lHini We also note that TBI have the 
"conven ient" property that both interactions of moderate 
strengtlPIEIl a,nd disorder'^'' '''* can sometimes enlarge the 
region of parameter space (for fixed spin-orbit coupling) 
where the topological phases exist, thus aiding their re- 
alization in experiment. 

In this paper, we focus on a particular two dimensional 
system: a tight-binding model on the square-octagon lat- 
tice illustrated in FiglT] that we show exhibits a number 
of Z2 topological phases. The square-octagon lattice con- 
sists of a square lattice with squares placed at the corner 
of each square plaquette in a square lattice. Part of the 
motivation to study this particular lattice comes from 
the knowledge that it supports topological phases in spin 
models, with the celebrated toric code^ as an effective 
low energy description. For certain parameter regimes 
in the spin models, Abelian and non-Abelia n p hases are 
obtained,S3 as well as a pseudo- Fermi surface.l^The con- 
nection between topological spin systems and topological 
band insulators was recently explored on the decorated 
honeycomb lattice in Ref.[3H]- We find that the square- 
octagon lattice tight-binding model of fermions we study 
indeed realizes a TBI phase over a wide range of param- 
eters, and there are a number of interesting quantum 
phase transitions appearing in the phase diagram. 

With this work, we expand the number of known 
lattices (honeycomb,^^ k ago me,^ checkerboard,ESl dec- 
orated honeycomb,'^ Lieb,'Hl pyrochlore,E2! perovskite,H3 
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FIG. 1. (color online) (a) (left) Schematic representation of 
square-octagon lattice with inter-square hopping t and second 
neighbor hopping A shown(dashed arrow line). Each vertex 
square and "half its emanating links is a unit cell of the lat- 
tice, (right) The Brillouin zone of the square-octagon lattice 
which has cubic symmetry (with a 4-point basis). High sym- 
metry points r, M, X and X' are labelled. (b,c) Illustration 
of gauge fields on nearest-neighbor hopping links. Different 
colors identify the way in which the gauge fields are applied; 
Blue (dark gray), green (light gray) and, red (gray) corre- 
spond to Uz,Ux and Uy gauge fields (as defined in the text), 
respectively. For each case a variety of inversion symmetric 
and asymmetric configurations are shown on the right. Ar- 
rows on each link indicate the sign of the gauge field. 



and diamondliESl) that support topological phases within 
a simple s-wave tight-binding approximation. When 
more orbitals are included, such as those with p and d 
character, the numb er of la ttices supporting topological 
phases is even larger .EEHSl Our main goal in this work is 
to help identify the conditions required to obtain topo- 
logical phases, and then predict how that physics can be 
observed in experiment. Our results are relevant to both 
solid state and cold atom systems. 

This paper is organized as follows. In SecllTlwe intro- 
duce a single-orbital model of fermions on the square- 
octagon lattice. We show topological phases are real- 
ized upon the inclusion of a second-neighbor spin-orbit 
cou pling term sim ilar to that introduced by Kane and 
Mele^. In Sec. Illll we show that in the absence of 
second-neighbor hopping, topological phases can be re- 
alized with non-Abclian gauge fields placed on the links 
of the lattice. As a function of the strength and type of 
non-Abelian gauge fields on the links, we find an excep- 
tionally rich phase diagram and study it in some detail. 
Finally, in Sec. |IV| we study the stability of the topo- 
logical phases to various symmetry breaking terms using 



the entanglement spectrum. In SeclVlwe summarize the 
main results of this work. 



II. KANE-MELE TYPE MODEL 
HAMILTONIAN 

The lattice we consider in this paper is shown in 
FigjT] (a). We first study a Kane-Mele type model 
Hamiltonian,l2El 



H = Hc, 



Hso, 
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The first and second terms in Hq describe the hopping 
of fermions with spin a on square plaquettes and be- 
tween neighboring plaquettes, respectively. We have set 
the nearest-neighbor on-plaquette hopping amplitude to 
unity and expressed all other energies in terms of this. 
Here c'^^ creates a fermion of spin a on site i and Cic an- 
nihilates a fermion of spin a on site i. The Hamiltonian 
Hso describes the spin-orbit coupling between second 
nearest sites whose relative position is encoded in the unit 



vector e ij, using the usual conventions: e*; 



|dl,x4|' 



the vector d\^ points from site j to a nearest neighbor 

site to both it and site ?, and d?- points from that near- 
est n eighbor t o the site i which is a second neighbor to 
site j.EEHnEil 

Exploiting the translational symmetry of the model, 
the Hamiltonian can be diagonalized as H = 

Ek.*L^k.*k., where *],, = (clk,,4k.. 4k.: 4kJ 
represents the four sites around the square unit cell (those 
that sit at the sites of the underlying square lattice) , and 
-ffko- is the Hamiltonian in Fourier space. The band struc- 
ture for different values of t and A of the model along 
the directions of high symmetry in the Brillouin zone is 
shown in FigMa-d). Each band is doubly degenerate 
due to the spin degree of freedom, and some bands touch 
each other at different crystal momenta in the Brillouin 
zone. For t ^ \ the behavior is shown in Figl2ja), and for 
A = we have a band touching at a single point, where 
a quadratic band touches a locally flat band, similar to 
what occurs in th e kagonie lattice and decorated honey- 
comb lattice .ISSESl An interesting feature of this lattice is 
that there are three bands (6 including spin degeneracy) 
that cross at the F and M points for i = 1 and A = (see 
FigJ2fb)), which are time reversal invariant momenta of 
the square lattice. At these points a locally flat band 
meets a Dirac-like structure coming from the other two 
bands. The situation is very similar to what occurs in 
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FIG. 2. (color online) Two top panels: the band struc- 
ture of the tight-binding square-octagon lattice model, Eq.([l|, 
along various high symmetry directions for the path shown in 
Fig[^a). Parameters are: (a) t = 0.8, A = 0, (b) t = 1, A = 0, 
(c) f = 1, A = 0.1 and (d) t = 1, A = 0.5. Bottom panel: phase 
diagram of the model at filling factor (e) 1/4 and (f) 1/2. 
Different phases are distinguished by colors as follows, blue 
(black): topological band insulator (TBI) state, yellow (light): 
band insulator (BI) and green (gray): semimetal (SM). Note 
that at 1/2 filing, there is no QSH state in this model and 
at 1/4 filling larger t tends to stabilize the topological phase. 
In (e) the solid vertical line A = is a semi-metal. The two 
intersecting dark lines in (f) running through the BI phases 
are lines where the system is a semi-metal with 4 Dirac points 
in the band structure as in Fig[2rd). 



the decorated honeycomb lattice model at hopping pa- 
rameters t' jt = 1.5 where a Dirac point is intersected by 
a flat band.^* 

In the present work we are interested in insulating 
phases, and one can follow different routes to gap the 
band structure of this model. One route is to localize 
the electrons to isolated square plaquettes. This is done 
by decreasing the hopping t in Hq. As i — >■ this model 
becomes gapped at 1/4 filling and 3/4 filling as shown 
in Fig[2]ja). It is trivial that at the extreme limit t = 0, 
the lattice becomes a set of disconnected plaquettes and 
is thus a trivial, non-topological insulator. All insulating 
phases created in this way are continuously (adiabati- 
cally) connected to this insulator. We therefore do not 
study this limit of Eq. ^ because this evidently gives rise 
to trivial insulators which are not the focus of this pa- 
per. The same trivial insulator phase also appears in the 
Kagome lattice by introducing a pattern of alternating 
bonds.SS 

Another approach to opening a gap is based on in- 
cluding a spin-orbit coupling term, Eq.(|3|, that changes 
the semi-metallic state of the model at 1/4 and 3/4 filling 
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FIG. 3. (color online) Energy dispersion of Eq.([l| in a strip 
geometry at different spin-orbit coupling A for fixed t = 0.5 
and 1/4 filling. Shown are the cases: (a) A = 0.2 (b) A = 0.5 
(c) A = 0.8. The states crossing the gap are edge modes, 
and on each edge of the strip right (solid red line) and left 
moving (dashed red line) states are degenerate at a single 
time-reverasal invariant momentum, which is a signature of a 
Z2 topological band insulator. Directly computing the Z2 in- 
variant indeed proves that (b) is a topological insulator state. 
In (c) there at also edge modes present, but the bulk is not 
ins ulatin g. The situation is reminiscent of that found in pure 
Sb.lSEZljji (d) we are considering a Haldane type model break- 
ing the time reversal symmetry by adapting a spinless model 
using the same parameters as in (b). The red (dark) and 
green (grey) lines stand for edge modes at right and left sides 
of the strip, respectively. The absolute value of the Chern 
number is 1 at 1/4 filling and 2 at 1/2 filling. 



into an insulator as shown in FigJ2[c). From a low-energy 
analysis around the Dirac -fflat band crossing points, one 
can deduce that a gap of value A — \AX + 2t — 2\ (note for 
i = 1 this vanishes for A = 0) opens up in the presence 
of spin-orbit coupling. However, at 1/2 filling, though 
the second and third bands separate from one another, 
for a wide range of spin-orbit coupling the Fermi energy 
crosses the bands resulting in a semimetallic state. In 
particular, at couphng |A/t| = 0.5, a Dirac crossing oc- 
curs at crystal momenta (±-|,±f) [see FigMd)], and 



upon further increasing of the coupling |A| a gap devel- 
ops. We will further discuss the implication of these Dirac 
nodes below. 

The gapped bulk phases at different filling fractions 
may be topologically trivial or nontrivial, and this is most 
easily seen via the edge modes.l^ the former is charac- 
terized by even number of Kramers pairs, while the latter 
is distinguished by an odd number of Kramers pairs. Be- 
cause our model has inversion symmetry, the Z2 invari- 
ants can be easily read off from the parity eigenvalues of 
the wavefunction at different time reversal momenta.'^Zl 
FigJ2je) and (f ) depicts how competition between param- 
eters in the Hamiltonian drives the model into different 
phases at 1/4 and 1/2 filling. The case of 3/4 filling re- 
sults in the same phase diagram as 1/4 filling because 
of the particle-hole symmetry present in the model (with 



A = 0). 

As discussed above, the topological insulator supports 
helical edge states at its boundary. To calculate the edge 
spectrum, we have considered our model on a strip ge- 
ometry which allows us to see the edge-projected mo- 
menta which clearly reveal edge states crossing the gap. 
Fig. [3] shows the spectrum at 1/4 filling for different val- 
ues of spin-orbit coupling corresponding to the different 
phases in Fig 2je) . Fig. Is] shows the behavior of edge 
modes through the transition between different phases. 
In the trivial insulator phase there are no gapless edge 
modes, effectively merged with the gapped bulk modes 
[see FiglsFa)]. However, as expected in the topological in- 
sulator phase [see FigjSJJb)], the edge modes traverse the 
bulk band gap. Importantly, there is an odd number of 
edge modes, which is the source of topologically robust- 
ness of these edge modes to weak disorder. Focusing at 
one of the edges of the strip, we can see that two travers- 
ing edge modes (solid and dashed red Hnes in Figjsj^b)) 
are a Kramers pair crossing each other at time reversal 
invariant momentum tt. The final figure of this sequence. 
Fig. Isle), is in a metallic phase, which is clear since the 
system is gapless. Some values of the edge-projected mo- 
menta have bulk conduction and valence bands merged 
into one another. However, even in this case we can see 
that there are traversing modes around the time-reversal 
invariant point tt. This mode can be shown to be local- 
ized on the edge of the system and behaves in many ways 
like the edge modes of the topological insulator phase. 
However, disorder on the edge can mix this state with 
bulk states, so it is not "topologically protected" .'^^'^ 
The physics of such boundary states in the presence of 
disorder has been investigated in Ref. HS] . In closing this 
part, we note that the stability of the topological phases 
to Rashba coupling and charge density wave order im- 
posed via a staggered on-site potential is very similar 
to that found for the honeycomb lattice^ and decorated 
honeycomb lattice.'^ We will investigate these effects in 
more detail in Sec |IV| via the entanglement spectrum. 

As we discussed earlier, the trivial insulting phase can 
not be adiabatically connected to a topological insulat- 
ing phase; along the way the gap A = must be closed. 
In this model, the low energy physics close to the gap 
closing point can be described by Dirac fermions. The 
corresponding Dirac node appear at the M point. At the 
phase transition the sign of the mass gap changes .1^02113 
Although this latter point often is a common feature of 
topological insulators and Hall systems which, respec- 
tively, preserve and break time reversal symmetry, one 
may wonder about the presence of only a single Dirac 
node. According to the Dirac low energy theory, each 
Dirac node should contribute half of the quantized Hall 
current, namely e^ /2h sgnfmV to the total Hall current 
at the edge of the system,E3 where sgn(m) is the sign 
of the mass gap. Thus, a naive expectation results in a 
half-odd integer when there are an odd number of Dirac 
fermions in the system, which contradicts the integer 
quantization of the Hall conductance as given by TKNN 



integer.^ In fact, the total contribution to the Hall con- 
ductance is not given by only summing up the individual 
contribution of Dirac fermions. Indeed, the TKNN in- 
teger (or Chern number) carries information regarding 
the Bloch states through the entire Brilloun zone, not 
just around the Dirac nodes. Thus, the vorticity con- 
tent of the Bloch states over the entire Brillouin zone is 
important.'"'*' However, the number of Dirac nodes gives 
the correct change of the Hall conductance through a gap 
closing process.ESl 

For a time reversal invariant system the Chern num- 
ber is zero. To get a nonzero value, one has to break 
time reversal symmetry. Here, we do this by considering 
only one spin species, say up, which is a kind of Haldane 
modepl^ on the square-octagon lattice. The bulk and 
edge modes at 1/4 and 1/2 filling for a strip geometry is 
shown in Fig. jsjd). The solid red (dark grey) and green 
(grey) lines traversing gap stand for edge modes localized 
at right and left of the strip, respectively. At 1/4 filling 
the Chern number in n = 1 which is consistent with the 
number of modes per edge's^ At this filling via the gap 
closing point the Chern number changes from n = to 
n = \, being equal to the number of Dirac nodes. 

As we have seen in Fig|2jf), at 1/2 filling there is no 
topological insulating phase. We can see that two Z2 
trivial insulator phases are separated by a gap closing 
point along two intersecting dark lines. At each gap clos- 
ing point, four Dirac nodes appear in the Brilloun zone 
[see FigJ2ld)]. For these trivial insulator phases, the Hal- 
dane type model considered shows the Chern number is 
not zero. Instead, it is n = 2 before and n = —2 after 
the gap closing point, which correspond to having two 
edge modes in the bulk gap as shown in Fig. [sld) (two 
red and two green lines). Thus, the Chern number, or al- 
ternatively the anomalous Hall current, will change by 4 
through the gap closing point. Once more we see that the 
change of the Chern number coincides with the number 
of Dirac nodes. 



III. HAMILTONIANS WITH NON-ABELIAN 

GAUGE FIELDS: TOPOLOGICAL INSULATORS 

AND TOPOLOGICAL PHASE TRANSITIONS 

In the preceding section we used second neighbor 
spin-orbit coupling to induce a topological insulator 
phase in a tight-binding model, which is a well-known 
paradigm for driving a transition to a topological in- 
sulator phase .f^I^EiElElHlI] Indeed, the essential ingredi- 
ent comes from the spin-dependent gauge fields^"^ (i.e., 
the spin-dependent second neighbor hopping in the most 
familiar cases). In this section we introduce a set of 
gauge fields living on the links between nearest-neighbor 
sites of the lattice. These gauge fields can be artificially 
inducecP^lt^ in the many-body Hamiltonian of ultracold 
atoms in optical lattices, where many rich behaviors can 
and haven been explored .123 For instance, ultracold atoms 
trapped in a honeycomb lattice, when subjected to non- 



Abelian gauge fields, exhibit various phases possessing 
different quantum orders^ due to the couphng of the 
gauge fields to the emergent relativistic quasiparticles. 
Here we address how such gauge fields may give rise to 
different phases of the model Eq.Q with only nearest- 
neighbor hopping. In particular, we are interested in 
whether such fields can stabilize a Z2 topological insula- 
tor phase. ^^ ®" Below we show that applying gauge fields 
to nearest-neighbor links leads to Zi topological band in- 
sulator phases as well as topological changes in the Fermi 
surface in the metallic phases of the model. 

We assume that the hopping terms in the first term of 
Eq.(^2^ are modified by some unitary matrices, t/, as,l^ 



H = 



E 

<i(y,j(j'> 



[U^j\„a,c\^C.ja' +^.C. 



(4) 



For spin 1/2 particles, a natural choice for unitary ma- 
trices will be the two-dimensional representation of the 
corresponding Lie group. 



U, 



^i-ia 



,U^ 



.Uy 



tis^y 



(5) 



where 7, a and /3 are parameters related to gauge fiuxes, 
and a"^ for (v = z, a;, y) stands for the usual Pauli ma- 
trices. Note that these gauge fields preserve the time re- 
versal symmetry, but may be applied in ways that could 
either preserve or break inversion symmetry. We consider 
two patterns for the modulation of hopping terms by ap- 
plying gauge fields on the links as shown in Figflrb,c). 
These two patterns give rise to rich phenomena that 
could potentially be realized in optical lattices. 



A. Gauge fields and topological insulating phases 

We start by focusing on the pattern shown in Figfljb). 
First we set a = /3 = implying that only hopping 
around square plaquettes on each point of the underlying 
square lattice have spin dependency. As depicted on the 
right hand of FigQb) , the gauge fields can be selected so 
that the inversion symmetry is preserved or not. Three 
different patterns are depicted. In the first pattern the in- 
version symmetry is broken. Here, the fluxes —27(27) are 
attached to square (octagon) plaquettes, and the degen- 
eracy between spin up and down is lifted. Only at time 
reversal invariant momenta do the different spin states 
remain degenerate. However, any nonzero value of 7 cre- 
ates a gap in the spectrum at 1/4 filling. This insulating 
phase is a topological band insulator as one pair of edge 
modes appear at the edge of the model in a strip geome- 
try, as shown in FigMa). In this case we can determine 
the spin-dependency of each edge mode, since the gauge 
field Uz conserves the spin up and down components. At 
each edge there in one Kramers pair crossing each other 
at time reversal momentum k^ = 0. In this figure the 
right (R) and left (L) edge are distinguished by red and 
green colors, respectively, and 'f (J,) stands for up (down) 
spin orientation. 



The second pattern of gauge fields in Figlljb) preserves 
the inversion symmetry of the lattice, and the square (oc- 
tagon) plaquettes carry —47(47) fiuxes. Different spin 
projections will have the same energy. Once again any 
nonzero value of 7 (except at 7 = 7r/4 where the flux 
pattern is equivalent to the zero-flux pattern because of 
the particle-hole symmetry) opens a gap at 1/4 filling. 
The main difference between the inversion symmetric and 
non-inversion symmetric flux patterns appears in nature 
of the edge modes. As shown in FigQb), there is one 
pair of edge modes traversing gap showing that for gauge 
fields that preserve the inversion symmetry the insulating 
phase is a topological phase. Direct evaluation of the Z2 
invariant also demonstrates this. When inversion sym- 
metry is broken, a topological phase still results, as seen 
in FigWa). However, the edge dispersion gets "split" 
relative to the inversion symmetric case. 

The third pattern of gauge fields in Figlljb) also pre- 
serves the inversion symmetry, but the plaquettes no 
longer carry fluxes. The model remains in the semimetal- 
lic phase at all fillings and for all values of 7, and behaves 
effectively as if no gauge fields are applied. This corre- 
sponds to the canceling part of the second-neighbor spin- 
orbit coupling on the square plaquetttes. In fact, accord- 
ing to Eq. (IT]) , the only nonzero contribution of spin-orbit 
coupling comes from the hopping between squares, while 
the second-neighbor hopping on the squares cancel each 
other. 

We also examined the stability of the topological insu- 
lator phase by turning on the gauge fields applied to the 
links connecting squares [red and green links in Figjljb)]. 
We found that as the value of either a or /3 is increased, 
there is a critical value for which the gap in the topolog- 
ical phase closes and therefore a quantum (topological) 
phase transition generally occurs. 



B. Topological phase transitions 

In this subsection we turn to a different method for 
modifying the hopping terms: applying the gauge fields 
as shown in FigfTJc). We assume that only the hopping 
terms around square plaquettes at the sites of the under- 
lying square lattice undergo such a modulation by the 
nontrivial gauge fields and we set 7 = 0. Here we can 
also consider different configurations respecting or vio- 
lating inversion symmetry. Let us first consider the case 
in which the inversion symmetry is preserved. Applying 
gauge fields with finite a and /3 gaps the spectrum at 1/4 
and 3/4 fillings. Again, we find that at 1/4 filling the 
gapped phase is topological. This fact is readily deduced 
by looking at the spectrum of the edge modes when the 
model is solved on a strip geometry, as seen in FigMc). 
The nonzero values of a and /3 breaks the symmetry be- 
tween up and down spins, and thus right or left movers 
are mixing states of different spin orientations. A direct 
evaluation of the Z2 invariant also shows this is a TBI. 
We note that at a = /3 = tt/2 the gap in the spectrum 











FIG. 4. (color online) Energy spectrum for the tight-binding 
model with nearest-neighbor gauge fields solved on the strip 
geometry. The gauge field modulations corresponding to the 
patterns shown in Fig[TJb,c) were used, (a) shows the edge 
modes with inversion symmetry broken and (b) inversion sym- 
metry preserved using the gauge fields in Figllrb) with param- 
eters a = /3 = and 7 = 7r/6 . (c) shows the corresponding 
edge modes with inversion preserved and (d) inversion bro- 
ken using the signs in FigfT^c) with parameters a = P = 7r/6 
and 7 = 0. The red (solid and dashed) and green (solid and 
dashed) lines stand for right and left of the strip, respectively, 
and arrows indicate different spin orientations. By counting 
the parity of the edge modes, we see that (a), (b) and (c) are 
topologically non-trivial, while (d) is topologically trivial. 



at 1/4 (and 3/4) filling vanishes and further increasing a 
and /? opens it again and returns the model to a topolog- 
ical insulating phase. 

A richer set of behaviors is observed when we consider 
the case where the gauge fields in Fig{l[c) break inversion 
symmetry. As we discussed before, there are two Dirac 
nodes (including spin) at 1/4 filling (when t ~ 1). Upon 
the inclusion of gauge fields, this degeneracy is split. To 
see this, let us consider the case of a = 0. Fig. [s] shows 
a set of band dispersions (left) alongside constant energy 
contours showing the evolution of the Dirac nodes within 
the Brillouin zone (right). This set of gauge fields cor- 
responds to the left side of the phase diagram shown in 
Figjs] (cyan dashed line Li), and the stars are different 
values of /3 in Fig. [5j Note that we have considered 
the band dispersions along the path connecting F to M. 
In the absence of gauge fields there is single Dirac node 
(doubly degenerate) at the center of the Brillouin zone, 
as seen in FigMa) and (b). As the gauge field /? increases 
from zero, two initially degenerate Dirac nodes are split 
apart and start to move along the path connecting center 
of the Brillouin zone to the its corner (F to M). These 
Dirac nodes are also intersected locally by a flat band. In 
the language of Ref. [S5] , there is no quantum phase tran- 
sition along this path as the total number of the Dirac 
nodes remains unchanged. They are only split and move 
about in the Brillouin zone. At the end of the path with 
/? = TT [see FigjsFi)] the two Dirac nodes are again merged 
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FIG. 5. (color online) Band dispersions (left) along F to M 
and corresponding contour plot (right) over the whole Bril- 
louin zone (dashed square) for a = 7 = 0: in (a,b) /3 = 0, in 
(c,d) /3 = ^/4, in (e,f) /3 = ^2, in (g,h) /? = 3^/4, in (i,j) 
P = n. Note that Dirac points move across the Brillouin zone 
as P evolves. The Dirac points are always intersected by a 
locally flat band. 



into a single one at the corner of the Brillouin zone. We 
could consider any border of the phase digram in FigjSJ 
they will present the same evolution. The scenario can be 
quite different if the gauge fields are applied in another 
manner. 

We now turn on both gauge fields a and (3 with a — jS 
[line L2 in Figis], still corresponding to the inversion 
asymmetric configuration in Figjljc) (because we already 
showed the inversion symmetric case results in a Z2 TBI 
when the spectrum is gapped). As seen in FigMa) and 
(b), the value a — (3 = OAw leads to a global gap at 
1/4 filling. The insulating phase here is a trivial insulat- 



ing phase since, as shown in FigllFd) , there are no states 
crossing the gap. However, further increasing of a = /3 
closes the gap at a — (3 — tt/4 [see Fig|6tc) and (d)] and a 
quadratic band touching point (QBTPpSj emerges along 
the line connecting F to M. In fact, in the parameter 
space spanned by a and /3, as shown in Fig[8J there are 
4 such QBTP indicated by solid circles. Since the gap is 
closed, we expect a quantum phase transition occurs at 
this point.133 The nature of this transition is character- 
ized by fine-tuning the gauge fields around the transition 
point. Fig. Wte) and (f ) reveal how bands cross each other 
(as a — /3 IS increase from the value of 7r/4) giving rise to 
a pair of Dirac nodes in the low energy dispersion along F 
to M. Thus, 8 Dirac nodes appear (in the low-energy sec- 
tor at 1/4 filling) in the Brillouin zone. Although part 
of the physics of the creation of Dirac nodes is similar 
to the non-Abelian hexagonal lattice,!^ the Dirac nodes 
here are created from a global vacuum. Once the nodes 
are created, they move away from each other towards the 
high symmetry points of the lattice, which are the F and 
M points. At a — (3 = 7r/2 four Dirac nodes at F and 
M merge into a single one. Thus, at this point we have 
two Dirac nods in the Brillouin zone as indicated by O in 
FigjS] As the gauge fields are increased further a converse 
process occurs. That is, the Dirac nodes at these points 
are split into four nodes and eventually at a = /3 = 37r/4 
are annihilated into vacuum. 

The two special cases of broken inversion symmetry de- 
scribed in Figjsland FigJGJfor the flux pattern in FigllTc) 
can be understood further. Motivated by the presence of 
two Dirac nodes in Fig. [5] and eight Dirac nodes in Fig. 
|6] (see also FigJ8|, it would be interesting to look for a 
path in the parameter space (a, /3) which connect those 
phases (indexed by the number of Dirac nodes) through a 
phase transition. We choose this path as /3 = — a-|-37r/4, 
and look at the deformation of the Fermi surface as the 
parameters are changed. This path is shown in FigjH] as 
L3, and for different values of gauge fields (triangles on 
L3) the Fermi surface is shown in Fig. [tI At a = the 
situation is already shown in Figjslg) and (h), where a 
Dirac like structure appears at the Fermi surface (for 1/4 
filling). At 1/4 filling, the Fermi level is illustrated by a 
red dashed line. As the gauge fields change, bands will 
cross the Fermi level filling and emptying fermi pockets. 
At very small values of a a finite Fermi surface (of holes) 
around the Dirac point appears, and becomes larger as 
the gauge fields is further increased. We see that a Fermi 
pocket (of electrons) is also formed in the Brillouin zone 
(between M' and F) . The structure of the Fermi surface 
will be deformed through a band inversion close to the 
Fermi level. This deformation leads to shrinking of the 
extended Fermi surface (around the Dirac point) to point 
like structures as shown in FigJTTc). It is interesting to 
note that through this evolution eight Dirac nodes are 
created at the Fermi level corresponding to FigJ6]jg) and 
(h). Thus, we see how a phase transition occurs between 
two phases each characterizing by different numbers of 
Dirac nodes. As the gauge field a increases further, the 
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FIG. 6. (color online) Band dispersions (left) along F to M 
and corresponding contour plot (right) over the whole Bril- 
louin zone(dashed square) for 7 = 0: in (a,b) a = P — O.lvr, 
in (c,d) a — P — 7r/4, in (e,f) a = P = O.Stt, in (g,h) 
a — P = 3n/8, in (i,j) a — P = 7r/2. Note the quadratic 
band touching point in (c) splits into 2 Dirac points in (e), 
(g), and (i). 



Fermi surfaces at different parts of the Brillouin zone are 
created and eventually at a = 37r/4, Fig. ^i), devel- 
ops two Dirac points as seen in Fig. [5](g) and (h). (The 
points a = 37r/4, (3 = and a = 0, /3 = 37r/4 have similar 
properties.) 

To summarize this section, our main results are the fol- 
lowing. We have studied the Hamiltonian (l2| with gauge 
fields applied to the nearest neighbor links, as shown in 
Eq. (41). Thus, in this section we have taken t = 1 m 
Eq.([2| and then apphed the gauge fields ^ to the finks 
using the patterns illustrated in Fig{l| (b j and (c) . We 
found that when the gauge fields preserve the lattice in- 
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FIG. 7. (color online) Band dispersions along the high sym- 
metry points of the Brillouin zone at difTerent gauge val- 
ues for inversion symmetry breaking flux patterns with /3 = 
—a + 37r/4: in (a) a = vr/8, in (b) a = 7r/4, in (c) a = Stt/S, 
in (d) a = 7r/2, in (e) q = 57r/8, and in (f) a = 37r/4. Along 
the evolution of a fermi pockets appear and disappear. 



version symmetry (they preserve time-reversal symmetry 
by construction) and a gap opens, the resulting state is 
a Zi TBI. Examples of TBI are shown in Fig. [I] (a) 
(which actually has inversion symmetry broken) , (b j and 
(c). The TBI obtained in this section can be compared 
with the TBI obtained in Sec|IT] where the gauge fields 
([5]) were absent and instead a second neighbor spin-orbit 
coupling (pi) was used to open a gap and drive the state 
into the TBI phase. Thus, if complex nearest-neighbor 
hopping parameters (that preserve time-reversal symme- 
try) are present a Z^ TBI can be stabilized. We speculate 
that it is possible to adiabatically continue a number of 
these states into each other (provided they have the same 
number of Dirac nodes on the edge) via a deformation of 
the Hamiltonian similar in spirit to that used in Ref . |38j . 

In the final part of this section we turned our attention 
to the flux patterns in Fig{l](c) that break inversion sym- 
metry. Here we found gapped states [some topological as 
shown in Fig. |4[a)] and gapless states depending on their 
strength. In the case of broken inversion symmetry we 
also found phase transitions in which the number of Dirac 
points in the Brillouin zone could change from 2 to 8, and 
we described these transitions in detail. We summarize 
them in the Fig. |8] Different parts of the parameter space 
can be distinguished by looking at the deformation of the 
Fermi surface. The border and center of the parameter 
space O represents the phase of the system with 2 Dirac 
nodes, and the thick blue lines represent the case with 8 
Dirac nodes. However, most region of parameter space 
is characterized by having a finite fermi surface. In fact, 
each band insulator phase will end up with a QBTP, and 
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FIG. 8. (color online) The parameter space of the lattice 
model with broken inversion symmetry gauge fields in Figjl] 
(c) for 1/4 filling. The yellow (light) thick lines indicate the 
trivial band insulator (BI) phases. The cyan dashed line (Li), 
thick yellow-blue line (L2) and black dashed line (L3) are 
paths we consider in FigJ5J Fig[6]and FigJTJ respectively. The 
thick blue lines correspond to phases with eight Dirac nodes, 
and solid circles indicate the quadratic band touching points 
(QBTP). The point O has 2 Dirac points in the Brillouin zone. 



these points (solid circles) evolve into eight Dirac nodes 
along the blue lines. 



IV. ENTANGLEMENT ANALYSIS OF THE 
TOPOLOGICAL INSULATING PHASES 

In this section, we further study the TBI phases via 
the entanglement spectrum.EHSil 'pjjg application of this 
metho d to t opological band insulators has been explored 
earlier jlSSEfil and we closely follow those works here. A 
key result of this paper is that a tight-binding model on 
the square-octagon lattice with only a single orbital per 
site supports TBI phases at 1/4 and 3/4 filling if second 
neighbor spin-orbit hopping of the Kane-Mele type is in- 
cluded (see Sec. Hi), or even if only first neighbor hopping 
is allowed provided non-Abelian gauge fields are added to 
the links of the lattice (see Sec. Ill I. We also noted that 



since the TBI obtained in each case has a single Dirac 
node on each edge, one expects them to be adiabatically 
connected via a deformation of the Hamiltonian similar 
to that used in Ref. [38 . 

In this section we return to the topological insulator 
phases discussed in Sec. [H] (similar results wouldbe ob- 
tained if we considered the models used in Sec. HI) and 



further characterize them by looking at the entanglement 



spectrum, 
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i.e. the spectrum of the reduced density 



matrix for a portion of the system. Let p be the pure 
density matrix of the whole system (obtained from its 
wavefunction) . Given p, we divide the system into two 
spatially distinct regions, A and B. The entanglement 
entropy is defined as Sa = — Tr[p^ln/9A], where pA is 
the reduced density matrix defined by tracing out the 
degrees of freedom of part B, that is, pA = TrB[/3]. 

The reduced density matrix pA is a versatile tool in 
both quan tum i nformation theory and condensed mat- 
ter physics I^^Ml xhe entanglement entropy Sa obtained 
from it is a useful measure of the quantum correlations 
between two parts of the system. Importantly, it ex- 
hibits different behavior on and off criticality.^^ ^° While 
it develops a logarithmic scaling in critical systems, it 
gets saturated for gapped syste ms an d exhibits the so- 
called area law behavior there.'^^'^ Since the ground 
states of topologically ordered states are highly entan- 
gled, the entanglement entropy also provides some infor- 
mation ab out th e topological order, namely the quantum 
dimension.ESESl 

However, the spectrum of the reduced density matrix, 
from which the entanglement entropy is extracted, con- 
tains richer information about topological phases-'^^™!! 
In this section we study the effect of various symmetry 
breaking perturbatio ns on the entanglement spectrum of 
topological insulators. ^^ ^^ This approach is complemen- 
tary to the study of the stability of topological band in- 
sulators discussed in Refs. |8I38I39| . 

A key difference between a topological insulator and a 
trivial insulator is tha t the former has "protected" gap- 
less edge modespESl t his fe ature will show up in the 
entanglement spectrum.ESESl Since the Hamiltonian in 
Eq.(ll| is noninteracting, the reduced density matrix of 
any part of the system can be f ully described in terms of 
correlation functions as follows,ESE3 
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where the matrix h is related to the correlation matrix 
G with elements Gij — {clcj) through the relation h — 
— log G+log(l— G). The partition function is determined 
by conditioning that Tr(p^) = 1 giving Z = l/det(l — 
G). The entanglement spectrum is then given by the set 
of eigenvalues of the matrix h. Thus, both h and G can 
be diagonalized in the same basis, and their eigenvalues 
are monotonically related to one another as^ 



o -& = o tanh( — ), 



(7) 



where gi and £/ are eigenvalues of correlation matrix G 
and matrix h, respectively. With this identification for 
the entanglement spectrum, we need only work with the 
spectrum of the correlation matrix. 

To calculate the entanglement spectrum, we consider 
a cut on the system that preserves the translational sym- 
metry along the cut. In this case, the momentum along 
the cut is a good quantum number and we can use it 
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FIG. 9. (color online) Illustration of the entanglement spec- 
trum for several different cases of Eq. ^ with t and A placing 
it in the Z2 topological phase: (a) both time reversed and 
inversion are preserved, (b) time reversal symmetry is broken 
by applying a Zeeman field while inversion is preserved, (c) 
time reversal is preserved and inversion is broken by Rashba 
coupling, (d) both time reversal and inversion symmetries 
are broken, (e) trivial insulating phase subjected to inversion 
symmetry presevering sublattice potentials, and (f ) trivial in- 
sulating phase with Rashba spin-orbit coupling. 



for labeling the spectrum. Thus, the single particle en- 
tanglement eigenvalues ei can be labeled by the linear 
momentum k along the cut, Si{k). In the actual calcu- 
lation, we consider two parallel cuts separated by many 
unit cells, and take the spectrum along one of the cuts. 
Note that this cut is not a physical cut, but only a cut 
separating regions where degrees of freedom are traced 
over to obtain the reduced density matrix pA- In fact, as 
shown in Refs. |65l55] the entanglement spectrum can be 
reconstructed from a spectrally flattened Hamiltonian. 

To study the stability via the entanglement spectrum 
of the Z2 TBI found in this paper, we consider the 
model Hamiltonian in Eq. (IT]) with the following perturb- 
ing terms: 



and 
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cLCo^a/3 X d^J)^Cjfi + H.c, (9) 



where the Hg and Hp^ describe the charge density mod- 
ulation and Rashba spin-orbit coupling. The 5i = ±<5 
can be chosen to either preserve or break the inversion 
symmetry of the lattice, while the Rashba coupling neces- 
sarily breaks the inversion symmetry. We consider both 
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cases to study the stability of the entanglement edge 
modes. As it turns out in the honeycomb lattice modeF' 
or its variant pS these perturbations suppress the topo- 
logical insulator phase at critical values in the present 
square-octagon lattice as well. 

First, we study the entanglement modes of the bare 
model without above perturbations. Fig. ^a) depicts 
how the edge modes of the topological insulator phase is 
reflected in the entanglement spectrum: the spectrum is 
gapless, just like the physical edge spectrum. The gapless 
physical edge modes are protected by time reversal sym- 
metry. However, the gapless modes in the entanglement 
spectrum are protected even if the time reversal symme- 
try is broken.ES To check this in the current model, we 
perturb Hamiltonian by a Zeeman term as 



H. = h^J2^l''o 



(10) 



which explicitly breaks the time reversal symmetry and 
gaps out the physical edge modes. However, as shown in 
Fig. l9[b) the gapless entanglement edge modes remain 
intact. This implies that the entanglement modes en- 
joy a higher degree of robustness than the physical edge 
modes. In fact, the entanglement edge modes are pro- 
tected by inversion symmetry r^^ For inversion symmet- 
ric systems, the inversion along the cut maps the right- 
hand part to the left-hand side of the cut. This invokes 
a kind of particle-hole symmetry,®^ and results in sin- 
gle particle entanglement energies having the property 
£r(~^) — —£i{k)- This relation clearly reveals that the 
edge modes of the entanglement spectrum at time rever- 
sal invariant momenta are degenerate with zero energy. 

However, the gapless nature of the entanglement edge 
spectrum can survive even if the inversion is broken, 
say by Rashba term in Eq.Q, provided the time re- 
versal symmetry is not also broken, as can be seen in 
Fig. ^c). Indeed, gapless physical edge modes lead to 
the de gener acy of the eigenvalues of the reduced density 
matrijpSMl and thus to the gapless entanglement edge 
modes. This follows from the argument that the 2^^ 
eigenvalues of the reduced density matrix Eq. (|6| can be 
written as ni=i[l + s;(e~^' — 1)] with s; = 0, 1. Trivially, 
a zero entanglement spectrum eigenvalue makes the den- 
sity matrix spectrum degenerate. On the other hand, ei 
can be related to the edge modes of a spectrally flattened 
Hamiltonian respecti ng th e ground state and phase of the 
original Hamiltonian.'SSESl j^ such a transformation, zero 
energy modes, which are often confined to the boundary, 
are reflected in the degeneracies of the eigenvalues of the 
reduced density matrix. 

Finally, breaking both time reversal and inversion sym- 
metry gaps out the edge modes of the entanglement spec- 
trum as shown in Fig. ^d). We see that the number 
of gapless edge modes in the entanglement spectrum is 
also consistent with the Z2 characterization of the topo- 
logical insulators.'^ To check the consistency with non- 



topological insulating phases, we consider two sets of pa- 
rameters corresponding to the trivial insulating phases. 
FigMe) and FigMf) illustrate the entanglement modes 
for trivial insulating phase resulting from the symmet- 
ric sublattice potential S — 1.2 and Rashba coupling 
A = 0.6, respectively (with t = 1) in Eq.(fl]). Both cases 
show an even number of edge modes consistent with their 
trivial Z2 values. 



V. SUMMARY AND CONCLUSIONS 

In this paper we introduced a single-orbital tight- 
binding model defined on the square-octagon lattice. 
First we studied a Kane-Mele type limit with spin-orbit 
coupling for second neighbor hopping, and showed that 
the model could support topological insulating phases 
depending on filling factors and the parameters of the 
model. In particular, we found that upon the changing of 
the spin-orbit coupling, the topological insulating phase 
turns into a semimetallic phase that could potentially 
have counterpropagating edge modes, which is qualita- 
tively similar to bulk Sb.l^S Then we considered a spe- 
cific modification of hopping terms in the Hamiltonian 
by coupling them to nontrivial gauge fields on nearest 
neighbor links (with no second-neighbor hopping). We 
showed such gauge fields can give rise to Z2 topological 
insulating phases. We found in all cases with inversion 
symmetry and a finite gap, the model is in a topological 
insulating phase, and in some cases with inversion sym- 
metry broken. Changing the gauge fields can also lead to 
topological phase transitions between different metallic 
phases. These metallic phases are characterized by the 
evolution of the Fermi surface and the number of Dirac 
nodes in the Brillouin zone. We also examined the sta- 
bility of the topological insulating phases by using the 
entanglement spectrum. While the time reversal break- 
ing perturbations gap out the physical edge modes, as 
long as the inversion symmetry is preserved, the entan- 
glement edge modes remains gapless. 

By providing a new example of a lattice that supports 
Z2 topological band insulators (realized via both second- 
neighbor spin-orbit coupling and nearest-neighbor non- 
Abelian gauge fields) in a simple s-band model, we have 
expanded number of known systems where these topo- 
logical phases exist. Our two approaches to achieving 
a topological state imply that both solid-state and cold 
atom systems can likely be found where the physics dis- 
cussed here is relevant. 
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